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Abstract 



We solve the Einstein equations in the Dvah-Gabadadze-Porrati 
model with a static, spherically symmetric matter distribution on the 
physical brane and obtain an exact expression for the gravitational field 
outside the source to the first order in the gravitational coupling. Al- 
though, this expression when confined on the physical brane reproduces 
the correct 4 — D Newtonian potential for distances ^ r ^ A, where 
A is a characteristic length scale of the model, it does not coincide with 
the standard linearized form of the 4 ~ D Schwarzschild metric. The 
solution reproduces the 5 — D Schwarzschild metric in the linearized 
approximation for r 3> A. 
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Linearized Analysis of the Dvali-Gahadadze-Porrati model 



Recently there have been several proposals which attempt to achieve localization of 
gravity. The Randall-Sundrum model [1] is based on the assumption that ordinary matter 
and its gauge interactions are confined within a four-dimensional hypcrsurface embedded in 
a five-dimensional space which is asymptotically anti-de Sitter. Since gravity is formulated 
by five-dimensional general relativity, the gravitational field produced by a mass source 
on the brane might deviate from the standard Schwarzschild metric. Calculation of the 
form of the metric produced by a static, spherical mass distribution on the brane, to first 
order [2] and to second order [3], [4], in the gravitational coupling revealed no tangible 
disagreement with the classical tests of general relativity. The Randall-Sundrum model 
was shown also to be consistent with the gravitomagnetic effect [5]. Various aspects of 
black hole physics in the framework of the brane world have been discussed in [6] . 



An alternative gravity localized model, the Dvali-Gabadadze-Porrati model [7] con- 
sists of a four dimensional hypersurfacc embedded in a five-dimensional space which is 
asymptotically Minkowski ( for a recent review of brane models see [8] . ) The main mo- 
tivation for this model is the competition between the bulk curvature scalar TZ and the 
corresponding intrinsic curvature scalar R on the brane. The action which describes this 
particular model is given by 

S = Ml J d'^xdyy/-g^^)n + M| j d'^xy/-g^^)R (1) 

where M5 stands for the 5 — D Planck constant and M4 is the 4 — D Planck constant. 
M5 and M4 are independent in general. The second term of (1) is included since matter 
localized on the brane might induce it through quantum fluctuations. This action yields 
the equations of motion 

MliTZmn - \ngran) + - ^-Rg^,)5!^5^J{y) = SXT^^^iy)- (2) 

In this paper we shall solve the equations (2), with T^^, ~ Mrjfj,i,S^{r) in the linearized 
approximation and we shall obtain an exact expression for the gravitational field outside 
the point source ( linearized analysis for the Randall-Sundrum model has been performed 
in [9].) For r 7^ 0, equation (2) implies 

n=~XR5{y), (3) 

and 

Tlmn + ^9mnRS{y) + X{R^, - ^Rg^,)SI^S:;,S{y) = (4) 



m2 



with A = being a characteristic length scale of the model. Here and throughout the 

5 

paper, we adopt the convention that the Greek indices take values 0-3 and the Latin indices 
0-4. The interesting feature of this scenario is that although the spectrum of linearized 
fluctuations does not include a massless graviton bound state which would correspond to a 
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4 — D graviton it reproduces the correct Newtonian interaction on the brane for a certain 
range of parameters [10]. 



In order to describe the real world, the Dvali-Gabadadze-Porrati scenario has to satisfy 
all the existing tests of General Relativity. The most general axially symmetric and static 
metric in Z) = 4 + 1 can always be brought to the following form 



ds^ = -e^dt^ + e^dr^ + e^r^d^^ + dy\ 



cJ2, 



(5) 



where dVt^ — dO^ + sin^ OdcjP' is the solid angle on S"^ and a, h and c are functions of r and 
y. The non-zero components of the A — D Ricci tensor for the general metric (5) in the 
linearized approximation are 
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while the components of the 5 — D Ricci tensor read 



(6) 
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(7) 



r 

where a' and d indicate differentiation with respect to r and y respectively. Consequently 
the A — D and the 5 — D Ricci scalars read 



R = a" + 2c" + -a' - -b' + -c' + ^(6 + c) 

rj^ rj^ rj^ rj^£i 

n = R+'d + b + 2'c. 



(8) 



These equations apply to the positive side of the brane, i.e., y > 0, the corresponding 
equations to the negative side of the brane, y < 0, are obtained by Z2 symmetry. The 
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Einstein equations (3) and (4) then demand that all of the components of TZmn vanish for 
y ^0 and dictate the following matching conditions at y = 

A 2 A 
a(r, y)|y=o = -^{a" + -a')\y=o + -^R\y=o 

b{r,y)\y=o = ~{a" + 2c" - -b' + -c%=o + ^R\y=o (9) 

c(r, y)|y=o = -^(c + -y + )ly=o + -^R\y=o- 

The yy component of the equation (4), together with equation (3) implies that i? = for 
y 7^ 0, which leads to i? = at y = by continuity. This simplifies the matching conditions 
(9) 

A 2 A 
a(r, y)|y=o = --{a" + -a')\y=o = -d\y=o 

kr, y)\y=o = -^(a" + 2c" - h' + ^c')|,=o = ^b\y=o (10) 

• / M , II 4,, a' -b\, A... 

c(r, y)|y=o = — 2~^^ + ~^ + —;—)\y=^ = 2^\y=o- 

The solution to the equation TZtt — reads 

1 f°° 

«(ny) = ^y dpp'^C{p)jo{pr)e~Py, y > 0, (11) 

where C{p) is a function of p to be determined and jo{x) is the spherical Bessel function. 
Substituting (11) into the first of equations (9) and taking into account that 6^{r) = 

I = ^ r dpp'joipr) we find 



Co 



As a result a(r, y) is given by the expression 

Co P^ 
2^^ Jo p+^Xp^ 

where Co is a constant to be determined by the requirement that the solution reproduces 
the correct Newtonian limit. The Newtonian limit is specified by the asymptotic behavior 
of goo on the brane for large r. Since p ^ ^ and p + ^p"^ ~ ;^ + 2^ we find that for r ^ A 



Co , , , Co 2GM , 
a(r,y = 0)^-^y^ dpjo(pr) = ^ = — (14) 



which determines the relation Co = — 47rGMA. 
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The yy-component of the equations of motion outside the brane provide the equation 
a + 6 + 2c = 0. If we impose the foUowing conditions to the solution 

hm a — hm h — hm c = (15) 

\y\,r^oo |y|,r^oo |y|,r^oo 

independent of the order of the hmits, we find that 

a + & + 2c=0, (16) 

which imphes i? = on the brane. By combining equations TZry = and (16) we derive 
the following expression for b{r,y) 

1 r 

Kr,y)^-:^ dr'r''^a{r' ,y). (17) 
Consquently by substituting (13) into (17) we find 

Kr,y) = / dp e y>0. (18) 

Ti" Jo p+^Xp^ pr 

Finally equation (16) provides an expression for c(r, y) 

c{r,y) = / dp \ JO -^^^ e P^ y > 0. 19 

Jo p+-^Xp^ pr 

Although a(r, y) satisfies the matching condition, 6(r, y) and c(r, y) fail to do so. The 
boundary condition for b{r, y) is (10) 

6(r,y)|^=o = -X{\a" + c" - ^6' + ^c')|,=o = ^Kr,y)|^=o. (20) 
By substituting equation (18) into (20) we find 

-Kr,y)|,=o - b{r,y)\y^o = ^^^^ = ^ ^^1) 

which is nonzero. The situation is reminiscent with the Randall-Sundrum model where 
two sets of solutions to the Einstein equations in two different coordinate systems, or 
equivalently, in two different gauges were constructed. The solution that we found is 
expressed in coordinates which are straight with respect to the horizon and is free from 
metric singularities far away from the source but it causes the brane to bend away from 
y = [11]. By transforming the solution to coordinates based on the brane we shall obtain 
another solution which satisfies the matching condition. In this set of coordinates the 
brane appears straight. We shall denote the solution in this gauge by a^{r,y), h^{r,y) 
and c^(r, y), where the superscript P means physical. 
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We shall now discuss the gauge symmetries of the Einstein equations, more specifically 
the coordinate transformations that respect the axially symmetric, static form of the five 
dimensional metric (5). Let's perform a coordinate transformation (gauge transformation) 
generated by v,u, functions of r and y such that r i— > r + v{r,y) and y ^ y + u{r,y) 
and demand that this particular coordinate transformation respects the form of the five 
dimensional metric (5). We find that the functions v,u obey to linear order the following 
relations 

v = 0, v' + u^0, (22) 

while the components of the metric transform under these residual coordinate transforma- 
tions as follows 

u 

Sa^O, Sb^2u', Sc^2-. (23) 

r 

We can easily verify that the Einstein equations remain invariant under the transformations 
(22). Thus, the physical solution will be given by 

= a + 5a, = b + Sb, = c + Sc. (24) 

We impose the boundary condition on on the brane, 

^b''{r,y)\y=o-b''{r,y)\y=o = 0. (25) 



2 

Under a coordinate transformation generated by u, v we find that 

\u (r, 0) - 2u (r, 0) = . (26) 

The parameters of the transformation, u and v, satisfy equations (22). Consequently, 

v = v{r), u{r,y) = -v'{r)y + x{r), (27) 
where x(r) is an arbitrary function of r. By substituting (27) into (26) we find that 

GMX , , GMX 

v{r) = u{r,y) = — -^y + x{r)- (28) 

The physical solution to the first order in the gravitational coupling takes then the form 
a [r,y) = / dp , o Jo(pr)e ^'^ 

Tl" Jo P+2^P 

,p, , 2GMX , p2 ) QM^^ 

Pf. .A - ^ r.^„r..^ _ ^1.-.. GMX^^ , 2x(r) 



c (r,y) = / dp A Mpr) ]e - -^r^y + 

Ti" Jo p+^Xp^ pr 2r-^ 
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Finally for 2GM <^r,y <^ \ the solution takes the form 



, , 2G'M, y + ir 
(^\Ti y) = : In 



4GM 



b{r, y) 



iirr 
GM 



y — ir 

y" 



Trr 



-tan - = 



2?/ 



2GM 2 1 r, 

(1 tan"^ -) 

r TT y 



(l-%)(l--tan-^-) + ^ 



c{r,y) = --{a + b) 



GM 



y 



Trr 



(30) 



2r 



n , y\f^ 2, _i r 2y 
(1 + ^)(1 tan -) 



TT 



y 



Trr J 



We observe that if we confine the metric on the brane y — 0, the physical solution reads 
for 2GM < r < A 



a^(r) 



2GM 



GM 



+ 2x', c^(r) 



GM ^ 2x 
2r r ' 



(31) 



where x is an arbitrary function of r. In order to recover the standard static, spherically 
symmetric metric on the brane, i.e. = at y — we choose x — This leaves 

b^ = , which differs from the standard form of the Schwarzschild metric and thus leads 
to different predictions for the bending of light. In the Randall-Sundrum model the choice 
of X that sets = leads also to b^ = which is the form of the Schwarzschild metric 
in the linearized approximation. 



Similarly for r,y ^ X 

, , 2GM\ 
a{r,y) ? 



J.2 _|_ 



, 2GMA, y ^ r, 
Kr,y) = ^:^{l-^tan-'-) (32) 

c(r-,y) = --(a + 6) = ^{ - -tan -). 

2 Trr^ y^ + r y 

We note that if we confine the metric on the physical brane y = 0, the solution reads 



„PM = _52^, t^(r) = ?^, c-M^O, (33) 



which agrees with the standard linearized form of the 5 — D Schwarzschild metric. 



In this paragraph we shall recapitulate what we have done in this paper. We solved 
the linearized Einstein equations in the model proposed by Dvali, Gabadadzc and Porrati 
with a static, spherically symmetric matter distribution on the physical brane and found an 
exact expression for the gravitational field. Although when confined on the physical brane 
the metric reproduces the 4 — D Newtonian potential for a certain range of parameters, 
more specifically for 2GM <^ r <^ A it fails to reproduce the standard Schwarzschild 
metric in its linearized form. Does this result then imply that the model is incompatible 
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with the experimental predictions of general relativity? Recently in [12] it was argued 
that the linearized approximation ceases to be reliable for this range of parameters and 
subsequently an exact classical solution to all orders in the gravitational coupling is needed. 
This argument is based on an early work by Vainshtein [13] in which he suggested that the 
limit of vanishing mass of massive general relativity may coincide with the predictions of 
massless general relativity in contrast with previous papers [14], [15] if the nonperturbative 
exact classical solution is used. The analog of the graviton mass on the brane is the inverse 
of the length scale A. The possibility then remains that a nonperturbative solution exists 
and it reproduces the standard Schwarzschild solution when confined on the physical brane. 
On the contrary the metric for r ^ A behaves as a 5 — D Schwarzschild metric. (See also 
[16] where it is claimed that four-dimensional behavior is recovered at large distances in 
the Dvali-Gabadadze-Porrati model.) We hope to report in the near future results which 
deal with these issues. 

Beyond the linear approximation, the properties of an event horizon which we estab- 
lished in a rigorous manner for the Randall-Sundrum model [17] apply to the DGP brane 
model as well. If a Schwarzschild black hole localized on the brane exists, the extension of 
the horizon into the bulk with the Gauss normal form of the metric (5) is tubular. This is 
not in contradiction with the analysis based on the linearized approximation, employed in 
this work, since the approximation breaks down for Ay ^ r^, thus permiting the horizon 
to extend to y y-^ oc. 
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